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ABSTRACT 

We analyze the X-ray spectra of the neutron-star (NS) X-ray transients Aql X-1 and 4U 1608-52, 
obtained with RXTE during more than twenty outbursts. Our aim is to properly decompose the 
spectral components and to study their evolution across the hard and soft X-ray states. We test 
commonly used spectral models and evaluate their performance against desirability criteria, including 
Lx oc evolution for the multicolor disk (MCD) component, and similarity to black holes (BHs) for 
correlated timing/spectral behavior. None of the classical models for thermal emission plus Conip- 
tonization perform well in the soft state. Instead, we devise a hybrid model: for the hard state a 
single-temperature blackbody (BB) plus a broken power-law (BPL) and for the soft state two thermal 
components (MCD and BB) plus a constrained BPL. This model produces Lx oc tracks for both 
the MCD and BB, and it aligns the spectral/ timing correlations of these NSs with the properties of 
accreting BHs. The visible BB emission area is very small 1/16 of the NS surface), but it remains 
roughly constant over a wide range of Lx that spans both the hard and soft states. We discuss impli- 
cations of a small and constant boundary layer in terms of the presence of an innermost stable circular 
orbit that lies outside the NS. Finally, if the BB luminosity tracks the overall accretion rate, then we 
find that the Comptonization in the hard state has surprisingly high radiative efficiency, compared 
to MCD emission in the soft state. Alternatively, if we assume that the radiative efficiency of a jet 
in the hard state must be less than the MCD efficiency in the soft state, while relaxing presumptions 
about the accretion rate, then our results may suggest substantial mass outflow in the jet. 
Subject headings: accretion, accretion disks — stars: neutron — X-rays: binaries — X-rays: bursts 
— X-ray: stars 



1. INTRODUCTION 

In low-mass X-ray binaries (LMXBs), a neutron star 
(NS) or stellar-mass black hole (BH) accretes mat- 
ter from a Roche-lobe filling, low-mass companion star 
through an accretion disk. X-rays are produced by the 
inner accretion disk and/or the boundary layer formed by 
impact of the accretion flow with the NS surface. The lu- 
minous and weakly magnetized NS LMXBs are classifled 
into atoll and Z source s based on their X-ray spectra l 
and timing properties (jHasinger &: van der KlislTl989f ). 
In a color-color diagram, Z sources trace out roughly Z- 
shaped tracks within hours to a day or so. Their X-ray 
spectra are normally soft in all three branches of the 
"Z", i.e., most of the flux is emitted below 20 keV. Atoll 
sources, however, show more dramatic spectral changes, 
albeit on longer time scales (days to weeks); their spectra 
are usually soft at high luminosities and hard when they 
are faint. 

Recently, it was found that some atoll sources can 
also exhibit Z-shaped tracks in the color-color diagram 
when they are observed over a large range of luminosity 
(|Muno et al.ll2002tlGierfinski fc Donell2002aD. How ever 
it was noted by several authors (iBa rret & Olive' '2OOI: 
van Straaten et all 120031 : iReig et all 2004; van der Klli 
2006[) that the properties of atoll sources (e.g., rapid X- 
ray variability, the order in which branches are traced 
out) are very different from those of the Z sources. The 
spectral states of atoll sources in the upper, diagonal 
and lower branches of these Z-shaped tracks are often 
referred to as the "extreme island" , "island" , and "ba- 
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nana" states/branches, respectively. However, in this pa- 
per we will use the terms "hard" , "transitional" , and 
"soft" states, respectively. 

The spectral modeling of accr eting NSs has been con- 
troversial for a long time (see iBarretl l200lL for a re- 
view). In the soft state, the spectra are generally de- 
scribed by models that include a soft/thermal and a 
hard/Comptonized component. Based on the choice of 
the thermal and Comptonized components, there are two 
classical models, often referred to as the Eastern model 
(af ter |M itsuda et al. J^989) and the Western model (af- 
ter |Wyte_et_alJll983). In the Eastern model, the ther- 
mal and Comptonized components are described by a 
multicolor disk blackbody (MCD) and a weakly Comp- 
tonized blackbody, respectively. In the Western model, 
the thermal component is a single-temperature black- 
body (BB) from the boundary layer and there is Comp- 
tonized emission from the disk. The color temperature 
for the thermal component (i.e., the temperature at the 
inner disk radius /cTmcd for a disk model or fcTbb for a 
boundary lay er model) is typica l ly in the range of ~ 0.5- 
2.0 keV (e.g..lBarret et a l. 2000; Oost erbroek et al.ll2001l : 
iDi Salvo et al.ll2000al:llaria et al. 2005[ ). For the Comp- 
tonized component, the same authors reported a plasma 
temperature of ^ 2-3 keV and a large optical depth of 
~ 5-15 in the soft state (for a spherical geometry). 

In the hard state, the spectra are dominated by a 
hard/Comptonized component, but a soft /the rmal com- 
pone n t is generally still required ^C hristia n fc Swank 
11997'; 'Barret et al."2000"; 'Chur ch fc Balucihska-Church 
HOOI; Gicrfinski & Done 2002b). The thermal compo- 
nent can be either a BB or a MCD, but the latter 
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seems to be ruled out by the inferred inner disk radii 
that are unphysically smaU. The color te mperature of 
this component is typically < 1 keV (e.g ., iBarret et alj 
[2OOI iChurch fc Balucihska-Churehl[200l . The inferred 
plasma electron temperature of the Comptonized compo- 
nent is typically a few tens of keV and its optical depth 
is ^ 2-3 in the hard state (for a spherical geometry). 
The hard state of atoll sources is considered by s ome au- 
thors to be associated with a steady jet (e.g., iFendeH 
I2006t iMigliari fc Fendeil I2OO6I ). The inverse Compton 
spectrum could then arise from the base of the jet, and 
synchrotron emission would contribute seed photons and 

erhaps a secondary contribution to the X-ray spectrum 

Markoff et'an[2005l ). 

In the past, various approaches have been made to 
further our understanding of X-ray spectra of accret- 
ing NSs. These include: (1) spectral surveys of a large 
number of sources, covering a wide range of luminosities 
(IChurch fc Balucihska-Churchll2001l : iChristian fc SwankI 
|1997[ ). (2) detailed studies of a large number of obser- 
vations from a single outburst of a NS X-ray transient 
jpicrhnski & Done 2002b; Ma ccarone fc Coppil l2003bt 
[Maitra fc Bailvn.2004). (3) Fourier frequency resolved X - 



rav spectroscopv (|Gilfanov et al.ll2003l : lOlive et al.ll2003D , 

and (4) comparisons of spectral and timing properties to 
those of BH LMXBs, to understand which features might 
be the result of the presence/ absence of a solid surfac e 
(|Wiinandsll200ll: lBarretl[200lt iDone fc Gierhhskil 120031) . 
However, a general consensus on the appropriate X-ray 
spectral model for the various subtypes and states of ac- 
creting NSs has not been achieved. 

In this paper we present an extensive study of a large 
number of observations of two NS transients Aql X-1 
and 4U 1608-52, using data o btained with the Rossi X- 
ray Timing Explorer {RXTE . iBradt et al.|[T993[) . during 
more than twenty individual outbursts. There are sev- 
eral advantages of using this archive. First, it allows us to 
compare the evolution of spectral properties of different 
outbursts systematically. Second, we can examine the 
behavior of a specific spectral component with changes 
in accretion rate, capitalizing on the large range in lu- 
minosity exhibited by these atoll-type transients. Third, 
compared with surveys using just a few observations but 
many sources, we can reduce the problems due to our 
poor knowledge of the parameters of the sources (e.g., 
the distance, inclination, and absorption). 

The goal of this paper is to find a spectral model that 
can well describe the X-ray spectra of NS X-ray tran- 
sients. Although the NS X-ray transients often change 
luminosity by several orders of magnitude and show sub- 
stantial diversity in outbursts, their well organized color- 
color and color-intensity diagrams clearly show that their 
spectral evolution tracks are narrow and thus repeatable. 
This compels the efforts to unlock spectral models for NS 
systems to explain the well-behaved sources in terms of 
accretion physics, and also to further investigate the dif- 
ferences and similarities between BHs and NSs. 

Aql X -1 and 4U 1608-52 have been classified as ato ll 
sources (iHasinger fc van der Klislll989tlReig et al.ll2000D. 
Recen tly^ iReig et aTT (|2004D and Ivan Straaten et al l 
(|2003l ) analyzed the general timing properties of these 
two sources. Here we will concentrate on their spectral 
properties. We describe our data reduction scheme in fj2] 
and show the long-term light curves and color-color dia- 



grams in Sj3l We perform and evaluate detailed spectral 
modeling in ^ In SJSJ we compare timing properties with 
BHs in order to further evaluate the spectral models. We 
argue that a particular model is most suitable for these 
NS transients, and then wc further explore the ramifi- 
cations of this model for X-ray states and the physical 
properties of accretions in ^Sl and 33 Finally we give our 
summary and discussion. 

2. OBSERVATIONS AND DATA REDUCTION 

For our analysis, we used all of the available RXTE ob- 
servations of Aql X-1 and 4U 1608-52 prior to 2006 
January 1. Data were ana lyzed from the Pr opor- 
tional Counter Array (PCA; iJahoda et al.l Il996f ) and 
the High Energy X-ra y Timing Experiment (HEXTE; 
iRothschild et al.l I1998D instruments. We utilized the 
best-calibrated detector units of each instrument, which 
are Proportional Counter Unit 2 (PCU 2) for the PCA 
and Cluster A for the HEXTE, and we extracted the 
average pulse-height spectra, one for each RXTE obser- 
vation. All spectral extractions and analyses utilized the 
FTOOLS software package version 6.0.4. Some standard 
criteria were used to filter the data: data of 20 seconds 
before and 20 seconds after type I X-ray bursts were 
excluded (see iRemillard et al.ll2006l ): the earth-limb el- 
evation angle was required to be larger than 10°; the 
spacecraft pointing offset was required to be < 0.02°. For 
faint observations, we additionally excluded data within 
30 minutes of the peak of South Atlantic Anomaly pas- 
sage or with large trapped electron contamination. 

We only considered observations that yielded data 
from the PCA, HEXTE and relevant spacecraft teleme- 
try. Only observations with PCA intensity (background 
subtracted) larger than 10 counts/s/PCU were used. We 
required the exposure of the spectra to be larger than 
five minutes for faint observations (source intensity lower 
than 40 counts/s/PCU) and two minutes for bright ob- 
servations. Appropriate faint/bright background mod- 
els were used when the source had intensity lower or 
higher than 40 counts/s/PCU. For the PCA, the spectra 
were extracted from "standard 2" data collection mode 
and the response files were created so that they were 
never offset from the time of each observation by more 
than 20 days. For the HEXTE, the program HXTL- 
CURV was used for spectral extraction, background 
subtraction, and deadtime correction. Finally, we ap- 
plied systematic errors of 0.8% for PCA channel 0-39 
about below 18 keV) and 2% for PCA chann el 40-128 
Krevkenbohm et al.ll200l iJahoda et al.ll2006f ). No sys- 
tematic errors were applied for HEXTE data. A sum- 
mary of the observations and other source properties is 
given in Table [H All further analyses and spectral fits 
were uniformly applied to each selected observation. 

3. LIGHT CURVES AND COLOR-COLOR DIAGRAMS 

The long-term light and color curves of Aql X-1 and 
4U 1608-52 are shown in Figure [T] This figure com- 
bines d ata from the PCA a nd the RXTE All-Sky Monitor 
(ASM: iLevine et anil996D . Typically, outbursts start in 
the hard state, evolve to the soft state, and finally return 
to the hard state during the decay. However, outbursts 
can be quite different from each other, e.g., in amplitude 
and duration. Moreover, in some outbursts the source 
did not enter the soft state. 
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Fig. 1. — Long-term light and color curves of Aql X-1 and 4U 1608-52, showing a large variety of outburst properties. Grey solid lines 
are from RXTE ASM and discrete symbols are from PCA observations representing different spectral states: hard (square), transitional 
(triangle), and soft (cross). The dashed line (hard color = 0.5) is a reference to help distinguish the different states. 
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Fig. 2. — Normalized color-color and color-intensity diagrams of 
Aql X-1 and 4U 1608-52 from ten years of pointed RXTE PCA 
observations. The PCA intensity is the normalized count rate from 
PCU2, and in these units the Crab Nebula yields 2520 c/s. The 
hard, transitional, and soft states are represented by blue squares, 
green triangles, and red crosses, respectively. The uncertainties are 
at 1 (T confidence and are normally smaller than the symbol size. 



Fig. 3. — The HEXTE intensity in two hard energy bands versus 
the PCA intensity at 3—20 keV. The solid lines show linear fits for 
these measurements. The hard and soft states appear increasingly 
different when they are compared using more widely spaced energy 
bands. 



For each observa tion we calculated X-ray colors, in a 
manner similar to iMuno et al.l (120021 ). Soft and hard 
colors were defined as the ratios of the background- 
subtracted counts in the (3.6-5.0)/(2.2-3.6) keV bands 
and the (8.6-18.0)/(5.0-8.6) keV bands, respectively. We 
normalized the raw count rates from each PCU with the 
help of observations of the Crab Nebula. For each PCA 
gain epoch, we computed linear fits (vs. time) to normal- 
ize the Crab count rates to target values of 550, 550, 850, 
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and 570 counts/s/PCU in these four energy bands. The 
normalized color-color an d color-interi s ity d iagrams in 
Figure [D resemble those in iMuno et al.l (|2002| ) , although 
our figure includes more data and we used a different 
normalization scheme for the softest energy bands. Al- 
though these two sources are atoll sources, their tracks 
in the color-color diagrams are Z-shaped, owing to their 
large range of luminosities. We combined the hard color 
and PCA intensity to pragmatically define the source 
states as shown in Figure [21 The hard state has a hard 
color > 0.65 (Aql X-1) or > 0.6 (4U 1608-52); the soft 
state has a PCA intensity > 400 counts/s/PCU and a 
hard color < 0.5 or a PCA intensity < 400 counts/s/PCU 
and a hard color < 0.36 (Aql X-1) or < 0.41 (4U 1608- 
52). All of the remaining observations are referred to 
as the transitional state. For all figures in this paper, 
the hard, transitional, and soft states are represented 
by blue squares, green triangles, and red crosses, respec- 
tively. We note that strong hysteresis is observed in both 
sources, that is, the hard-soft transition generally occurs 
at hi gher X-ray flux com pared to the soft-hard transition 
(e.g. jMaccarone fc Cop pi 2003a). The transitional-state 
observations are mostly from the decay phases of the out- 
bursts, since the rise is often sparsely covered. 

Figure [3] shows the relation of source intensities in more 
widely spaced energy bands. The top panels are the 
HEXTE 20-50 keV intensity versus the PCA 3-20 keV 
intensity. The most striking aspect of these panels is that 
there are two nearly linear tracks corresponding to the 
soft and hard states. The bottom panels are the HEXTE 
50-150 keV intensity versus the PCA 3-20 keV intensity. 
We can still see the linear track in the hard state, but 
the sources are generally not detected above 50 keV in 
the soft state. 

While the light curves in Figure [T] show substantial 
diversity in outburst amplitude and duration for a given 
source, Figures [2] and [3] show that the superposition of all 
observations on color-color and color-intensity diagrams 
shows well organized spectral states and very common 
behaviors in these two sources. The aim of this paper is 
to capitalize on this organization and the comprehensive 
RXTE data archive, in order to determine the best way 
to model the X-ray spectra across these different states. 

4. SPECTRAL MODELING 

4.1. Spectral Models and Assumptions 

In this work, we fitted the X-ray spectra of Aql X- 
1 and 4U 1608-52 with several different models. The 
PCA and HEXTE pulse-height spectra were fitted jointly 
over the energy range 2.6-23.0 keV and 20.0-150.0 keV, 
respectively, allowing the normalization of the HEXTE 
spectrum, relative to the PCA spectrum, to float between 
0.7 and 1.3. For soft-state observations, HEXTE spectra 
were used up to 50 keV because the flux at higher photon 
energies was negligible (Figure [3|) . 

In the typical description of the NS continuum spectra, 
some form of pure thermal radiation is often combined 
with another radiation process, commonly presumed to 
be some form of Comptonization, although synchrotron 
radiation might also be involved (§1). Hereafter, we sim- 
ply refer to emission other than the pure thermal radi- 
ation as the "Comptonized" component. Two forms of 
thermal radiation were considered: BB and MCD mod- 
els (bbodyrad and diskbb in XSPEC respectively). The 



BB model provides the color temperature (fcTbb) and the 
apparent radius (i?bb; isotropic assumption) of the BB 
emission area, while the MCD model provides the appar- 
ent inner disk radius (i?mcd) and the color temperature 
at the inner disk radius (fcTmcd)- 

As for the modeling of the Comptonized component, 
we considered both a broken power-law model (bkn- 
power in XSP EC, hereafter BP L) and the Comptoniza- 
tion model bv iTitarchufl (fl99l (CompTT in XSPEC). 
The CompTT model computes the Comptonization of a 
Wien input spectrum of " seed photons" by a hot (single 
temperature) plasma with a uniform covering geometry. 
On the other hand, the BPL component can be consid- 
ered as a functional approximation for Comptonization 
under complex conditions or in combination with another 
radiation process like synchrotron radiation. We gave 
these two models an equal opportunity to handle the ef- 
fects of Comptonization when we tested different kinds 
of spectral decomposition. Thus BPL or CompTT was 
combined with BB and/or MCD in a variety of ways, 
as summarized in Table [2J The BPL model has four 
parameters: two photon indices, a break energy (i?b) 
and a normalization parameter. The CompTT model is 
parametrized by the seed-photon temperature (kTs), the 
plasma electron temperature (kTc), the optical depth (t), 
a parameter describing the geometry of the Comptoniz- 
ing cloud (either spherical or disk- like), and a normaliza- 
tion parameter. 

All models also included a Gaussian line. Its central 
line energy was constrained to be between 6.2-7.3 keV, 
targeting the Fe Une 6.4 keV (Asai ct al. 2000). The av- 
erage best-fitting value was ^ 6.6 keV for both sources. 
The intrinsic width of the Gaussian line (ct) was fixed 
at 0.1 keV. This is consistent with the ASCA result of 
IChurch fc Balucihska-Cluirehl (pOOll ). The PCA has en- 
ergy resolution ~ 1 keV, limiting the need for a precise 
value. An interstellar absorption component was also in- 
cluded with the hydrogen column fixed at A^h = 0-5 x 
10flcm"2 foj. Aql X-1 (jChurch fc Balucinska-Churcl] 
200lh and 1.0 x 10^^ cm'^ for 4U 1608-52 (JPrnminx et al.l 
Since no eclipses or absorption dips have been 
observed, the two sources are likely not high-inclination 
systems and we assumed their binary inclinations to be 
60° . We scaled the luminosity and radius related quanti- 
ties u sing distances of 5 kpc for A ql X-1 ([R utledge et alj 
I2001h and 3.6 kpc for 4U 1608-52 (jNatalucci ct al.,,.200Q) . 
unless indicated otherwise. 

4.2. Comptonized + Thermal Two-component Models 
4.2.1. The Problem of Model Degeneracy 

First we consider models with two continuum com- 
ponents: one is Comptonized and the other is thermal 
(Models 1-4 in Table [2]). It turns out that the spectral 
fitting is inherently non-unique. To illustrate the con- 
tribution of each component and the problem of model 
degeneracy, we use two observations of Aql X-1. One is a 
hard-state observation on 2004 February 21 (hard color 
0.73, intensity 178.5 counts/s/PCU), the other is a soft- 
state observation on 2000 October 27 (hard color 0.33, 
intensity 1262 counts/s/PCU). 

Figures [4] and [5] show the unfolded spectra of these 
two observations using different Comptonized + ther- 
mal models. All models give acceptable spectral fits, de- 
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Fig. 4. — The unfolded spectra and residuals of two sample 
observations of Aql X-1 using different kinds of models. One is in 
the hard state and the other in the soft state. The panels for models 
with CompTT are from spectral fits with seed-photon temperatures 
< 0.5 keV (cold-seed-photon models). The total model fit is shown 
as a solid line, the Comptonized component as a dotted line, the 
thermal component as a dashed line, and the Gaussian line as a 
dot-dashed line. 

spite the fact that the spectra each contain well over 10^ 
counts. Figure [4] shows that there is not only degener- 
acy from the choice of the thermal component (i.e., BB 
or MCD), but also from the choice of the Comptonized 
component (i.e., BPL or CompTT). Moreover, even if we 
choose Model CompTT-t-MCD or CompTT+BB, there 
are still two competing minima, one with best-fitting 
kTs < 0.5 keV and the other with best-fitting kT^ > 1 
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Fig. 5. — The unfolded spectra and residuals of two sample 
observations of Aql X— 1 with best-fitting kTs > 1 keV using Model 
CompTT-(-MCD (hot-seed-photon model). The total model fit is 
shown as a solid line, the Comptonized component as a dotted line, 
the thermal component as a dashed line, and the Gaussian line as 
a dot-dashed line. 

keV. Their corresponding unfolded spectra can be seen in 
Figures |4] and [5l respectively. Table [3] gives the detailed 
results using the CompTT-hMCD and CompTT+BB 
models for our two representative spectra. Hereafter, 
we call models with best-fitting fcTg ^ 1 keV "hot-seed- 
photon models" and models with best-fitting fcTg < 0.5 
keV "cold-seed-photon models". Hot-seed-photon mod- 
els require that the modeling of the Comptonized compo- 
nent takes into account the spectral curvature expected 
from having the seed photons close to the observed band- 
pass (Done et al. 200^). We note that the temperature 
of the chosen thermal component increases and its nor- 
malization decreases significantly when the seed-photon 
model flips from the hot to cold solution, for a given 
observation. 

Examination of the other observations shows that this 
seed-photon problem is quite general. CompTT-l-MCD 
and CompTT+BB models (Table [2]) yield two solu- 
tions each: photon temperature kTs ^ 1 keV and 
kTs ^ 0.5 keV. The inferred parameters of the ther- 
mal component are also quite different (Table [3] and 
references following). For the soft-state observations, 
the inferred temperature of the thermal component 
is typicall y > 1 keV for the cold-seed-photon mod- 
els (e.g., [ Barret ct^ I200GI : lOosterbroek et al.l 1200 It 
IWhite et a l. f988). Otherwise, if the hot-seed-photon 
models are used, the temperature of the thermal compo- 
nent is normally < 1 k eV and thus less than fcTg (e.g., 
iDi Salvo et aLll2000allbl : [laria et al.ll2005D . For the hard- 
state observations, if BB is used as the thermal compo- 
nent, the size of the BB emission area is typically very 
small, ^ 2 km, fo r the cold-seed-photon models (e.g. , 
Barret et al.l 120 031: 'Ch urch fc Balucihska-ChurchI 120011 : 
Gierhnski fc Do ne 2002^^ However, with the hot-seed- 
photon models, the BB emission a rea can be compa- 
rable to the size of the NS (e.g.. iBarret et al.l 120001 : 
iNatalucci et~aLl l2000t iGuainazzi et al.l 119981) . Replace- 
ment of BB with MCD in cold-seed-photon models yields 
very small inner disk radii as already found by other 
authors (see references in §1). We also point out that 
many of the references cited above made use of broad 
band spectra from BeppoSAX. This means that the seed- 
photon problem is also present for instruments that 
have extended low-energy spectral coverage (see also 
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iFarinelli et al.l (|2005D ). We also analyzed the BeppoSAX 
observations of Aql X-1 and 4U 1608-52 and found the 
same problem to be present. 

The low-energy limit of PCA is 2.3-3.0 keV, depend- 
ing on the gain setting epoch. We cannot constrain kTg 
in the cold-seed-photon models because the peak energy 
flux of the seed-photon spectrum (Wien approximation) 
is at 3kTs < 1.5 keV. On the other hand, for the hot- 
seed-photon models, the fits to PCA spectra produce a 
cool thermal component in the soft state making it is dif- 
ficult to constrain the temperature and normalization of 
this component, and sometimes the thermal co mponent 
is even not required (|GierHnski fc Don^l2002bl ). 

4.2.2. Efforts to Resolve Model Degeneracy 

As outlined in the previous section, there is a degener- 
acy in X-ray spectral models for NS LMXBs. The prob- 
lem is based on the fact that acceptable fits {xf, criteria) 
can be obtained for either the hard or soft states by using 
any combination of two ambiguous components: a ther- 
mal spectrum (BB or MCD) plus a Comptonized compo- 
nent (BPL, hot-seed CompTT, or cold-seed CompTT). 
The various models convey (very) different pictures for 
the structures and energetics of NS accretion. 

In our descriptions of various model details, we have 
begun to mention arguments that have been offered to 
evaluate competing models in terms of physical implica- 
tions derived from the fitted spectral results. This strat- 
egy of performance-based evaluations of spectral models 
is most effective when there is a self-consistency issue at 
stake. Multiple evaluation criteria will help us to choose 
a spectral model that may be superior in overall suit- 
ability. Useful considerations can involve single parame- 
ters, spectral evolution, and the relationship between the 
hard and soft states. The comparison between the atoll 
sources and BHs can provide further constraints on this 
choice. Below we list one such consideration used in the 
literature, followed by three considerations offered in this 
paper. We describe them in term of problems that are 
encountered in a particular model/spectral state, and we 
track these problems in the fourth and seventh columns 
of Table El 

1. "R" problem: the inner disk radius in models 
with MCD is too small, i.e., much less than the 
size of the NS (IChurch fc Balucinska-Churchl l2001t 
iGierlihski fc Donell2002b[ r 

2. "L" problem: £mcd oc T^^^ (i.e., constant in- 
ner disk radii) for MCD component is not sat- 
isfied for any meaningful range of luminosity in 
the soft state; We note that BH systems do show 
Lined cx T^^H when they are in the soft (thermal) 
state (|Kubota fc Donel[200l . 

3. "T" problem: kTg > fclbb in models with BB; if we 
attribute the BB component to the boundary layer 
and accept that the boundary layer should have a 
highe r temperature than the dis k (|Mitsuda et al.l 
llQSl iPopham fc SunvaevI [200lh . kT^ seems un- 
likely to be larger than fcTbb. 

4. "P" problem: the Comptonization fraction is not 
consistent with the power density spectrum, as- 
suming that a comparison of atoll sources and BHs 



in timing properties is relevant. We will explain 
this in more detail in Sj5l 

In the definition of the "L" problem, we use the phrase 
"over some range of luminosity" in recognition of the 
possibility that the disk may deviate from the MCD's 
geome tric assumptions, as do BH accretion disks at high 
L,„cd (jKubota fc Dond 12004( 1. To help evaluate the "R" 
problem, we inferred the NS radii of our two sources from 
the spectral fitting to Type I X-ray bursts (denoted as 
Rburst)- During the decay phase of some bursts, an ap- 
proximately constant burst emission area is derived when 
temporal series of the burst are fitted with a BB model 
after subtracting the average persistent emission to iso- 
late the burst from other radiation components. Such an 
emission area is expect ed to be roug hly similar to the ap- 
parent size of the NS (jLewin et al.lll993f ). Our fits using 
PCA data gave ii!burst ^ 8 km for Aql X-1 and ~ 7.2 km 
for 4U 1608-52 (at distances in Ta ble [J), which are con- 
siste nt with previous results (e.g., iKovama et al.l [l98ll : 
Nak amura et al.lll989l ). However, we note that these val- 
ues are given without any corrections, e.g., hardening 
factors and surface occultation by the inner disk. These 
issues are further discussed in SjH 

4.2.3. Model Fit Results 

We fitted the observations of Aql X-1 and 4U 1608-52 
using Models 1-4 in Table H The PCA does not al- 
low us to gain meaningful constraints on the properties 
of the (cool) thermal components using the hot-seed- 
photon models, as explained above. Thus, for models 
with CompTT, we just show results of cold-seed-photon 
solutions. In these models, fcTg kT^ {kT^ is typically 
2-3 keV in the soft state and several tens of keV in the 
hard state), and fcTg is far below the PCA energy limit. 
Thus, the Wien approximation of the input seed photon 
is valid. 

The top four rows in Figures [6l [7| and [8] show the fit 
results for the commonly used Comptonized -I- thermal 
models. They correspond to Models 1, 2coid, 3, and 4coid 
in Table [21 where we also list the mean xt values. The 
errors represent 90% confidence limits for a single pa- 
rameter, and the circled points are those with relatively 
large errors. There is not only strong similarity between 
these two sources, but also strong similarity between the 
results of these four models. The extreme difference of 
the hard and soft states can also be seen in these plots. 

Figure |6l shows the variation of the color temperature 
of the thermal components (fcTbb for BB and fcTmcd for 
MCD) with the PCA intensity. All four models show 
that the temperature of the thermal components of both 
sources increases with intensity in the hard state. In 
contrast, the temperature is almost constant (for models 
with BPL) or decreases (for models with CompTT) with 
increasing luminosity in the soft state. 

Figure [71 shows the luminosity of the thermal compo- 
nent versus its color temperature, fcTbb or fcT^cd- For 
reference, we also show the lines for constant radius, 
assuming Lx = AnR'^aT^. The NS radii (t gX^j) are 
shown with dotted lines in this figure. The dashed 
lines correspond to i? = 1.9 km and 1.3 km for Aql X- 
1 and 4U 1608-52, respectively; they are derived from 
the fit to the BB radius values obtained from Model 
6 (see below). With increasing luminosity in the re- 
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Fig. 6. — Variation of the temperature of the thermal component 
(feTbb for BB and fcT^cd for MCD) with the PCA intensity. Cold- 
seed-photon solutions are shown for Models with CompTT. 

spective thermal component, R 

mcd Rhh incrGcises in 
the soft state. For the case of MCD-related Models 1 
and 2coidi this behavior in the soft state contradicts the 
basic prediction of the accretion disk model and war- 
rants the "L" problem in Table [H Besides, from the 
comparison with the i?burst lines, we find that the in- 
ner disk radius i?mcd in the hard state and some part 
of the soft state is simply too small, even after taking 
into account expected correction factors (see j}6l) . This is 
con sistent with the results of iGie rlihski fc Pond ()2002bl ) 
and lChurch fc Balucihska-Chureh (|2001i ). Thus, the "R" 
problem also applies to these two models. 

For BB- related Models 3 and 4coid, the luminosity evo- 
lution of the effective radius is tied to the evolution of 
the boundary layer, which is much less certain. This 







1 


3 BPL + BB j 




.1 




































D ■■ □ 










1 

0.1 

1 

0.1 [ 

0-^ ■ 
0-^ ■ 



4 CompTT + BB 






- ^ 



































5 MCD+BB 
MCD 



® 





kT.cd or kT,, (keV) 

Fig. 7. — The luminosity of the thermal component versus its 
characteristic temperature. Models 5 and 6 produce luminosity 
evolution for BB and MCD components in the soft state (red 
crosses) that are nearly parallel to reference lines, which show 
Lx with constant emitting surface area. The dotted lines 

correspond to the NS burst radii ( m.'l.'ll and the dashed lines cor- 
respond to i? = 1.9 km and 1.3 km for Aql X-1 and 4U 1608-52, 
respectively, assuming Lx = ii^R^crT'^. 



is why the "L" problem is not applied to the BB com- 
ponent in Table [2l If we attribute the thermal compo- 
nent BB to the boundary layer, it would imply that the 
boundary layer is measured with almost constant surface 
area in the hard state, but it spreads out with constant 
color temperature in the soft state. One possible ex- 
planation for the constant color temperature during the 
spreading of the boundary layer is that the local flux 
reaches the Eddington limit. The spreading layer model 
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Fig. 8. — The luminosity of the Comptonized component ver- 
sus the luminosity of the thermal component(s). The dashed and 
dotted lines are for reference. They connect the points where the 
ratio of the luminosities of the thermal component and the Comp- 
tonized component is 1.0 and 0.07, respectively. The latter is a 
typical value for the hard state. 



imanov fc Poutanenl I2OO60 predicts that the color 
temperature of the boundary layer is ^ 2.4 keV for a NS 
with i?ns = 15 km and M = IAMq and solar composi- 
tion of the accreting matter (smaller i?ns would predict 
higher value). A problem with this interpretation for 
our results of Model 3 or 4coid is that the critical value of 
the temperature for area expansion is ~ 1.5 keV. Higher 
temperatures are only seen during Type I X-ray bursts 
from these two sources with maxima normally > 2.5 keV. 
Thus, the critical value of the temperature ~ 1.5 keV 
from Models 3 and 4coid does not match expectations for 
boundary layer spreading. We will show later that these 
two models are also disfavored from strong similarities of 
the timing properties between BH and NS LMXBs. 

In Figure [8] we investigate the luminosity evolution of 
the Comptonized component (BPL or CompTT) versus 
the thermal component (s) (BB or MCD or their sum if 
both are used) . For the thermal component the bolomet- 
ric luminosity is used while for the Comptonized compo- 
nent we integrated from 1 keV to 200 keV for CompTT 
and from 1.5 keV to 200 keV for BPL. The 1.5 keV lower 
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Fig. 9. — The unfolded spectra and residuals of two sample soft- 
state observations of Aql X-1 using Model 6. The total model fit is 
shown as a solid line, the MCD component as a dotted line, the BB 
component as a dashed line, the CBPL component as a dot-dot- 
dashed line, and the Gaussian line as a dot-dashed line. The high- 
Lx observation is the same as the sample soft-state observation 
seen in Figures |4] and [5] The low-Lx observations is on 1999 June 
3 (hard color 0.32, intensity 243 counts/s/PCU). 

limit for the BPL integration is chosen so that this com- 
ponent does not extend below the temperature of the 
MCD in the soft state, i.e., when the BPL is steep and 
the lower limit matters. The errors of the luminosity of 
BPL, CompTT and MCD primarily depend on the uncer- 
tainties in their respective normalizations. In Figure[5]we 
also show two reference lines: dashed and dotted. They 
connect the points where the ratio of the luminosities 
of the thermal component and the Comptonized compo- 
nent is 1.0 and 0.07, respectively. The latter is a typical 
value for the hard state. In both the hard and soft states, 
the luminosity of the Comptonized component increases 
with the luminosity of the thermal component, but the 
two states follow different tracks. In the soft state, the lu- 
minosities of these two components are relatively close to 
each other, while in the hard state, the luminosity of the 
thermal component is only ^ 10% of the Comptonized 
component. 

The uncertainties for the thermal-component luminos- 
ity from the hot-seed-photon models are quite large from 
PCA data in the soft state. However, hot-seed-photon 
models with BB as the thermal component have the "T" 
problem ( tj4.2.ip . Hot-seed-photon models with MCD as 
the th ermal component have the "L" problem in th e soft 
state ()Gierhhski fc Dond l2002bl : iDone et al.l l2002f ) and 
the "R" problem in the hard state from our investiga- 
tion (not shown). 

Since the luminosity evolution of the thermal compo- 
nent in the soft state is either in violation of the basic 
model (MCD) or, at best, suspicious (BB), we contin- 
ued to investigate alternative models such as those in 
the next section. 

4.3. Double Thermal Models 

Early analyses of NS spectra in the soft state con- 
sidered the possibility that we might detect thermal 
components from bo th the disk and boundary layer 
(|Mitsuda et all I1984D . This "double thermal" model 
(i.e., MCD-I-BB or Model 5 in Tabled was applied to 
the soft-state observations of Aql X-1 and 4U 1608-52, 
and the results are shown in the fifth row of panels in 
Figures [6] and [T] Data points for this model are omitted 
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Fig. 10. — The fraction of Comptonized luminosity versus the 
total luminosity, using Model 6. Ledd = 1-8 X 10"^** erg/s, assuming 
Mns = IAMq for both sources. 

when is > 2 (only for this model). It turns out that 
this model works very well in the soft state when the 
luminosity is high, but with the decreasing luminosity, 
x1 becomes large. Either source gives a mean xt ~ 1-2 
and 3.5 for observations with source intensity > 800 
counts/s/PCU and < 500 counts/s/PCU, respectively 
(also see Table [2]). However, Model 5 is remarkably suc- 
cessful in its implication that the disk and the boundary 
layer both remain at constant sizes for a substantial range 
of luminosity (Figure [7]). Moreover, the inferred size of 
the boundary layer is close to that inferred by Models 3 
and 4coid for the hard state. The failure of Model 5 for 
soft-state observations at lower luminosity is apparently 
due to small levels of Comptonization, since the fit resid- 
uals are pointing chronically positive at photon energies 
above 15 keV (PCA data). 

Given the interesting luminosity evolution for the soft 
state with Model 5 (no Comptonization), we added a 
weakly Comptonized component, while continuing to as- 
sume that both the MCD and BB are visible in the soft 
state. The difficulty here is how to model weak Comp- 
tonization by adding a third component. Freely adding 
a BPL or CompTT component is obviously not feasible, 
because BPL or CompTT plus one thermal component 
is sufficient to model the entire spectrum as shown by 
Models 1-4. As a first attempt, we tried several forms 
of constrained CompTT (e.g., couple kTs to kT^h or 
fcTincd; force kTc > 10 keV; or do both), but the fits 
always yielded large fractions of Comptonization for low- 
luminosity observations in the soft state. This problem 
could be due to RXTE's lack of low-energy coverage. 

On the other hand, we found that the entire soft state 
remains weakly Comptonized if we adopt the following 
constrained BPL (CBPL) model: the break energy is 
fixed at 20 keV (the best-fitting break energy is typ- 
ically > 20 keV in the hard state and is ~ 15 keV 
in the soft state from Models 1 and 3), and the ini- 
tial photon index is required to be < 2.5 (a typical ini- 
tial p hoton index in the soft state of BHs is also about 
2.5: iRemillard fc McClintockl I2OO60 . Therefore, we de- 
fine Model 6 as follows: the hard state is still modeled 
by BPL+BB (no BPL constraints; same as Model 3 for 
the hard state) and the soft/transitional states are mod- 
eled by MCD-FBB-hCBPL. Figure [HI shows the unfolded 
spectra of two soft-state observations using this model, 
one at low luminosity and the other at high luminosity. 

The bottom row of panels of Figures El E] and [5] show 
the results obtained for Model 6. Figure [5] shows that 
with the increase in the intensity, fcTbb increases both in 



the hard and soft states, with the tracks that are clearly 
separated. The BB temperature reaches a maximum ^ 
2.5 keV for Aql X-1 and ~ 3.0 keV for 4U 1608-52, simi- 
lar to the peak tem peratures of Type I X-ray bursts from 
these two sources jK ovama et al.|[l981l : iNakamura et al.l 
[1989: Galloway et al. 2006) . The temperature at the in- 
ner disk radius kT^cd increases from ~ 0.5 to 2.0 keV, 
also correlated with the intensity. 

In Figure [71 we see that Model 6 remarkably produces 
results where Lx oc T'* for both the MCD and BB com- 
ponents in the soft state. Furthermore, the inferred emis- 
sion areas of the boundary layer from the hard and soft 
states have essentially the same value. The emission area 
of the boundary layer is small compared with the size of 
the NS. We realize that the size of the BB emission area 
slightly increases with decreasing luminosity in the hard 
state, especially for 4U 1608-52. However, we note that 
the BB curve lies well below the line of i?buist at all lu- 
minosities. 

In Figure [51 we show the luminosity of the Comp- 
tonized component versus the total luminosity of the 
double thermal components. The behavior in the soft 
state for Model 6 is quite different from that in the 
first four models. Here, with increasing luminosity of 
the thermal components (i.e. Lbb + imcd)j the lumi- 
nosity of the Comptonized component first increases and 
then decreases. At the highest luminosities, the Comp- 
tonized component is negligible, as implied by the suc- 
cess of Model 5 in the same region. Figure [TOl shows the 
fraction of the Comptonized luminosity versus the total 
luminosity. It looks similar to the color-intensity dia- 
gram in Figure [2l suggesting that the hard color tracks 
the degree of Comptonization fairly well. 

It should be noted that there are other kinds of weak- 
Comptonization approximations like constrained power- 
law (photon index < 2.5) or constrained cutoff power-law 
(photon index < 2.5 and cutoff energy > 10 keV) that 
gave results for the soft state that are similar to those 
obtained using CBPL. For the hard state, Figures 6-8 
show that the Model A:coid (CompTT-|-BB) gives simi- 
lar results to Model 6 in terms of the properties of the 
thermal component and the fraction of Comptonization. 
Model 6 is successful, but there are no claims that it is 
either a unique solution to the problem, nor an adequate 
depiction of Comptonization other than the estimate for 
the fraction of energy related to Comptonization. 

5. TIMING PROPERTIES AND COMPARISON WITH 
BLACK HOLES 

Compared to Models 1-4, the use of a double thermal 
model -I- CBPL (i.e.. Model 6) for Aql X-1 and 4U 1608- 
52 produces dramatically different results for the lumi- 
nosity evolution of the thermal components in the soft 
state and for the implied significance of Comptonization 
in the soft state. 

There are strong similarities in timing pro perties be- 
tween the BH and NS LMXBs (see, e.g., IWiinandsl 
(l200l '). which can be used to further assess these dif- 
ferences. In Figure [TTl we show the integrated root- 
mean-square (rms) power in the power density spectrum 
(0.1-10 Hz and energy band 2-40 keV) versus the hard 
color for Aql X-1 and 4U 1608-52. The two sources show 
very similar timing properties. The rms, normalized to 
a fraction of the source's mean count rate, is very small 
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Fig. 11. — The integrated rms power in the power density spec- 
trum (0.1-10 Hz and energy band 2-40 keV) versus the hard color. 
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Fig. 12. — The integrated rms power versus the fraction of lu- 
minosity contained in the thermal (MCD) component for 6 black 
hole systems. 

(< 0.05) in the soft state, and the values increase with 
the hard color progressing through the transitional and 
hard states. The rms versus hard-color relations from 
Fi gure | 11 | can be compar e d with similar plots for BHs in 
iRemillard fc McClintockl ()2006D (panels g in their Fig- 
ures 4-9). Based on such a comparison we conclude 
that, when considered in a model-independent manner, 
the way in which the strength of the X-ray variability 
changes as a function of spectral hardness is very similar 
for BHs and our two NS transients. Since the hard color 
effectively traces the fractional contribution of the ther- 
mal and Comptonized components to the X-ray spectrum 
in BH systems, it is interesting to see which of our spec- 
tral models conserves this similarity when rms is plotted 
versus the fractional contribution of the thermal compo- 
nent (s). 

In Figure [12] we pres ent such a plot for the same 6 BH 
binaries analyzed by IRemillard &: McClintockl (|2006[ ). 
Two different integration limits were investigated for the 
thermal and Comptonized luminosities, but the overall 
behavior was insensitive to our exact choice. It can be 
seen that for BHs in the soft state, the luminosity trac- 
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Fig. 13. — the rms versus the fraction of luminosity contained 
in BB and/or MCD components as evaluated for different spectral 
models. 

tion of the thermal component is quite high (> 75%) 
and the rms is very small < 0.05. For Aql X-1 and 4U 
1608-52 we plot the rms versus the fraction of thermal 
luminosity for Models 1, 3, and 6 in Figure [T3l Note that 
in case of Model 6 the sum of MCD and BB was used 
to determine the fractional contribution of the thermal 
components. Comparing Figures [12] and [13] one can see 
that only Model 6 reproduces the same dependence of 
rms on fractional thermal luminosity for Aql X-1 and 4U 
1608-52, as for the six BHs. In Models 1 and 3, the lumi- 
nosity fraction of the thermal component can be as low as 
40% and 20%, respectively, while continuing to show very 
low values of rms power. Similar results are obtained for 
Models 2 and 4 (cold and hot included for both cases). 
Hence, for the models outlined in Table [2] only Model 
6 links weak Comptonization with low rms power, as is 
clearly evident in the properties of BHs. Models 1-4 for 
the soft state are therefore described as having the "P" 
problem (see g^iland Table [1]). 



PHYSICAL INTERPRETATIONS OF MODEL ( 
THERMAL COMPONENTS 



FOR 



Since Model 6 has many attractive advantages over the 
other models, we further explore its implications in terms 
of the physical properties of the accretion flow in both 
the hard and soft states. 

There is a well known difficulty in deriving true radii 
from the apparent dimensions (i.e., i?bb, ^burst, and 
^mcd extracted from model components of the X-ray 
spectrum BB, BBburst, and MCD, respectively). Nev- 
ertheless, the small value of i?bb and the nearly constant 
values of i?bb and i?nicd across a large range in luminosity 
motivate cautious efforts to discuss physical interpreta- 
tions. 
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If we assume the BB emission area to be a latitudinally 
symmetric equatorial belt, then the BB emission area of 
the belt should be 

Aboit = 4^iVbbi^?okpc/^fc(*, S) W, (1) 

where iVbb = ^bb/^iokpc normalization of the BB 

component (isotropic assumption), / represents a spec- 
tral hardening factor (sometimes expressed as the ratio 
of the color to effective temperature) , and fc is a geomet- 
rical correction factor taking into account the emitting 
geometry and any occultations by the accretion stream. 
The k factor depends on the disk inclination (i), the lat- 
itude range (from the NS equator) of surface emission 
(S). We note that this factor is also strongly dependent 
on the properties of the occulting accretion stream. Un- 
less otherwise indicated, we assume a geometrically thin 
but optically thick accretion stream so that the half belt 
on the other side of the accretion stream is invisible to 
the observer. To have some sense of this k factor, we 
give two examples: k{i = 60°, S = 10°) = 1.69, and 
k{i = 60°, 6 = 90°) = 1.33. From this perspective, the 
area of the NS, e.g., for Aql X-1 in this study is: 

A„3 = 4^iVbursti^?okpc/^fc(*, S - 90°) km^ (2) 

assuming that the entire NS surface is radiating (using 
the asymptotic value of -Rburst late in the burst). This 
expression conveys the difficulty in gaining accurate in- 
ferences of NS sizes from X-ray burst measurements. The 
raw values for -Rburst given in ij4. 2. 21 ignore / and k. 

Perhaps of greater interest is the effort to understand a 
key measurement result of this paper: -Rbb/-Rburst — 0.25 
(-^bb/-^burst — 1/16) for Aql X-1. If we assume that the 
/ values cancel in these different applications of the BB 
model and use the fact j4beit/^ns = sin 5, we obtain 

iVbbfc(i,<5) _ k{i,S) 



sin (5 = 



(3) 



iVburstfc(i,(5 = 90°) I6k{i,6 = 90°)' 

In this case, the Rhh/ Rburst value is equivalent to an 
equatorial belt with a half- angle of i5boit — 9° for i = 30°, 
or (5boit — 6° for i = 60°. We stress that this scaling 
estimate ignores the annular half width of the occulting 
stream (^stream) itself. A more realistic estimate is (5bcit ^ 
^stream + 6° for i = 60°, implying fc > 1. 

Regardless of the details of 5, we must confront the 
implication of measuring A'bb/A^burst ~ 1/16, noting that 
A'burst is the asymptotic value that should screen out ef- 
fect of the radius expansion and momentary disruption 
in the occulting stream f ij4.2.2p . How is it possible that 
accreting NSs could maintain an almost uniformly small 
area of BB emission through hard and soft states that 
span such a wide range in luminosity (i.e., 0.005 to 0.5 
-^edd)? Such simple results suggest a picture in which a 
geometrically thin accretion stream feeds a rather well- 
defined impact zone, where the accreting gas may effi- 
ciently release energy before spreading over the remain- 
der of the NS surface. Since the scale height of the inner 
accretion disk is expected to increase with the accretion 
rate, it is difficult to understand our results without the 
help of an inner-most stable circular orbit (ISCO). As 
illustrated for BHs, the effective potential of an ISCO 
creates a pinch on the vertical structure of the accretion 
stream that can s uppress variations in the scale height 
of the inner disk (|Abramowicz "eTaniigTSh . To the ex- 
tent that Model 6 remains viable with further scrutiny. 



the small size of _Rbb over a large range in luminosity 
should be examined as a means to infer that these NSs lie 
within their ISCOs. This may provide another observa- 
tional link to general relativity, while providing potential 
constraints on the NS equation of state. 

The ability of Model 6 to restore the expected lumi- 
nosity evolution of the MCD component motivates efforts 
to compare the MCD radius and relative luminosity to 
the results determined for the BB component. In Fig- 
ure [71 it is evident that the raw values of -Rmcd are only 
slightly less than the values of i?burst (dotted line) for 
Aql X-1 and 4U 1608-52. Given the uncertainty in the 
different correction factors that must be applied, respec- 
tively, to the disk and the burst radii in order to derive 
physical sizes, our results may still be consistent with ex- 
pectations framed by the preceding discussion, i.e., that 

Rns ^ R'mco ^ -Rdisk ■ 

The comparison of BB and MCD luminosities is a trick- 
ier topic. The apparent luminosity of the boundary layer 
(Lbb) is about one third of that of the disk (imcd) in 
the soft state, neglecting the weak Coniptonization. In 
contrast, we might expect Lbb ^ ^mcd if the boundary 
layer is inside the NS ISCO, as sketched above, where 
the accreting material may impact the surface with a 
large relative velocity that includes a component along 
the radial path. We note, however, that while imcd is es- 
timated in a true bolometric sense, this is not the case for 
Lbb- We need to correct Lbb by k{i,5) (Equation[T]) due 
to the geometry of the boundary layer and occultation by 
the accretion stream. For Aql X-1, it is ~ 1.7 (i = 60°) 
or ~ 2.6 [i = 30°), and these factions can become much 
larger if ^stream ~ <5bcit, as noted above. We conclude 
that there is considerable uncertainty in our final results 
as to whether the total energy losses at the boundary 
layer are less than the total bolometric luminosity of the 
other spectral components. 

7. PHYSICAL CONSEQUENCES OF MODEL 6 FOR THE 
HARD STATE 

Figure [14] shows the luminosity evolution with time 
during a well-covered outburst of Aql X-1 in 2000. The 
light curve of this outburst is shown on the top panel. 
The middle panel shows the luminosity evolution of each 
spectral component, using our hybrid Model 6, where 
the BPL and MCD dominate the apparent luminosity of 
the hard and soft states, respectively. The bolometric 
BB luminosity (Lbb) curve resembles the total luminos- 
ity curve, and it provides a reference measurement that 
allows us to further compare the hard and soft states. 
The ratio of BPL (isotropic) and MCD (bolometric) lu- 
minosities (Lbpi + Lined) to Lbb is shown in the bottom 
panel. We note that the values of this ratio are elevated, 
in part, by the k factor described for Eq. 1. Nevertheless 
the high luminosity of the hard state, relative to Lbb, is 
apparent in the bottom panel of Figure [TH 

To compare the hard and soft state luminosities in a 
more global manner. Figure [TSl shows (Lbpi + Lmcd) ver- 
sus Lbb for all of the RXTE observations of Aql X-1 and 
4U 1608-52. The behavior of these two sources is quite 
similar. The states are largely separated by vertical lines 
near (uncorrected) Lbb ~ 10'^^ erg/s. Here, the hard- 
state values for (Lbpi + Lmcd) are clearly higher than 
those for the soft state, by a factor of ^ 6. For the hard 
state, the vertical axis represents Comptonization, since 
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Fig. 14. — The luminosity evolution of different spectral com- 
ponents (middle panel) during the 2000 outburst of Aql X-1, as 
viewed with Model 6. Blue, green, and red symbols denote hard, 
transitional, and soft states, as done previously. Note that the 
MCD component (points connected by a solid red line) is not di- 
rectly visible with the PCA in the hard state. The top panel shows 
the PCA count rate (2-30 keV). The bottom panel shows the ra- 
tio: (Lbpl + ^mcd) / J^bht S'lid the value is clearly highest in the 
hard state, when Comptonization appears to dominate the X-ray 
spectrum. 
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Fig. 15. — The luminosity of non-BB components (i.e., 
BPL-I-MCD) versus the luminosity of the BB component. At the 
transitional BB luminosity that separates hard and soft states, the 
hard-state BPL appears to be very luminous compared to the soft- 
state MCD. The two reference lines, dashed and dotted lines, con- 
nect the points where the ratio of the luminosities of the BB com- 
ponent and the non-BB components is 1.0 and 0.07, respectively. 



90% of the total luminosity in the hard state is due to 
Lbpi- In the soft state, the vertical axis is effectively the 
disk luminosity. On the horizontal axis, we assume that 
Lbh represents emission from the boundary layer, and we 
can then use the nearly constant value of i?bb across the 
hard and soft states to motivate the presumption that 
Lhh quantifies the mass accretion rate onto the NS sur- 
face. Using these assumptions, there are still different 



ways to interpret Figure [T5l 

First, if we further assume that the total accretion rate 
(M) flows into a boundary layer with constant 6hcit, then 
Lbb effectively tracks (M), and we can compare hard and 
soft states in terms of their radiation efficiency. Unlike 
Lhh, imcd and Lbpi are not expected to require large 
correction factors, although the origin of the BPL is very 
uncertain. We would then conclude from Figure [15] that 
Comptonization in the hard state has higher radiative 
efficiency than thermal disk emission in the soft state, 
and this occurs across a wide range in the accretion rate. 
This idea is a radical departure from expectations. 

Alternatively, we can refrain from making assumptions 
about M, and choose a line of constant (Lbpi + imcd) in 
Figure 1151 Quite generally, then, one would conclude 
that less mass reaches the NS surface in the hard state, 
compared to the soft state, for a given level of luminosity 
on the vertical axis. This might suggest that the major- 
ity of the mass flow is disrupted from the normal path 
to the NS in the hard state. Given the expectation for 
radio je ts in the hard state of accreting NS LMXBs (e.g., 
iFender 2006; M igliari k Fend er 2006), this could further 
imply that a significant fraction of mass accretion in the 
hard state is redirected to outflow in the jet. The energy 
source for this scenario could be problematic. 

Finally, it is also possible that some of our preliminary 
assumptions are incorrect. Changes in accretion geome- 
try across the state boundary, e.g., from correlated jumps 
in ^bcit and 5 stream, would be one way to invalidate the 
assumption that Lbb uniformly tracks the NS surface ac- 
cretion rate across both states. All of these ideas need to 
be investigated for other sources and in different ways. 
8. SUMMARY AND DISCUSSION 

We have analyzed ten years of RXTE observations 
of two NS X-ray transients, Aql X-1 and 4U 1608-52, 
through more than twenty individual outbursts. Al- 
though there is much variety in these outbursts, the spec- 
tral and timing properties are well organized , as has been 
demonstrated in previous investigations (|Muno et al.l 
[2002; Gicrhnski & Done 2002a). It is also wefl known 
that several different models may adequately fit the X- 
ray spectra of NS LMXBs, while each one may convey 
very different physical implications, e.g., the tempera- 
tures and sizes of the emitting regions. We therefore 
attempted to design performance-based criteria for eval- 
uating these models, hoping to gain insights on a model's 
self consistency, while taking advantage of the large data 
archive for these two atoll-transient prototypes. 

Among these criteria, there is an important quantity, 
which is the energy fraction for Comptonization, i.e., how 
much energy is directly visible as pure thermal radiation 
and what remaining fraction is diverted and expressed as 
inverse Compton radiation. This quantity is interesting 
from an accretion-energetics point of view. It can be 
calculated from any of these models, and furthermore, 
it is the gateway for comparing NSs and BHs in terms 
of the relationship between energetics and the measured 
timing properties. 

First we examined the classical two-component models 
with a Comptonized component (BPL or CompTT) and 
a thermal component (MCD or BB). We confirmed that 
many combinations can describe the X-ray spectra suc- 
cessfully. Moreover there are two possible seed-photon 
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solutions for models that utilize CompTT, a problem 
that is also present for instruments with extended low- 
energy coverage like BeppoSAX. In the soft state, the 
two-component models behaved poorly for a variety of 
reasons shown in Table [21 Figure [71 and Figures [T2HT31 
We then progressed to consider a double-thermal model 
(MCD-I-BB) for the soft state, and we found a need to 
account for weak Comptonization at the low-Lx end of 
the soft-state track. We then devised a hybrid "Model 
6" in which the hard state is modeled by BB+BPL while 
the soft state is modeled by MCD+BB+CBPL, where 
the third component is a constrained BPL (i?b = 20 keV 
and Fi < 2.5) that can assume the role of weak Comp- 
tonization. 

It turns out that Model 6 offers three great advantages 
over the classical Comptonized + thermal Models. (1) It 
produces an Lx oc relation for both the MCD and BB 
in the soft state with a sufhciently large value for the in- 
ner disk radius (compared to the NS radius inferred from 
Type I bursts). (2) The emission area of the BB extends 
from the hard to the soft state with essentially the same 
value. (3) The fraction of Comptonization inferred by 
this model is also consistent with that for BHs for given 
values of rms power in the power density spectra. When 
BH binaries are in the thermal state, the Comptoniza- 
tion fractions in the energy spectra are very low and the 
rms power is also very low, rms < 0.05. Atoll sources in 
the soft state show similar rms values, but only Model 6 
implies low Comptonization for the soft state. 

With regard to previous investigations of NS X-ray 
spectra (§1), our Model 6 splices together the Western 
model for the hard state (BB -I- a heavily Comptonized 
disk) and a three-component model for the transitional 
and soft states (BB -1- MCD -I- constrained BPL) that 
represents a BB plus a very weakly Comptonized disk. 
The soft-state portion of this model was heavily influ- 
enced by the original version of the Eastern model (Mit- 
suda et al. 1984). We further note that low Comptoniza- 
tion solutions (r ~ 1; kTc > 20 keV) for the soft state 
were derived with a two-component Eastern model for 
4U 1608-52 with Tenma data (Mitsuda et al. 1989), and 
these are the closest results that we can find to those 
(soft state) of Model 6. 

The ad hoc scheme in Model 6 to handle weak Comp- 
tonization with the CBPL component for the soft state 
is obviously a topic that needs further work. Model 
6 for the soft state is a three-component continuum 
model. We did not try to test three-component mod- 
els in an open-handed manner. Instead, we focused on 
low-Comptonization solutions, progressing from the suc- 
cess of Model 5 (no Comptonization). The CBPL used 
in this paper is just one possible solution; our investi- 
gation showed that constrained power-law with photon 
index < 2.5 or constrained cutoff power-law with photon 
index < 2.5 and cutoff energy > 10 keV would also pro- 
duce similar results to those obtained with the CBPL. As 
other alternatives, we tried to add a constrained version 
of CompTT to the two thermal components, by coupling 



kTs to /cTbb or fcr,„cd, or requiring kTc > 10 keV, or do- 
ing both. However, the fits always yielded very strong 
Comptonization in the low-luminosity soft state, but the 
fit residuals for Model 5 suggested that only weak Comp- 
tonization is present in the low-luminosity soft state. 
This could be due to the limited low-energy coverage 
of RXTE data. We investigated BeppoSAX, which has 
extended low-energy coverage, for possible new insight. 
Unfortunately, each X-ray source has only one soft-state 
observation from BeppoSAX, and both of them are at 
high luminosity. Future observations with broad-band 
instruments like Suzaku could be important to this is- 
sue. As for the hard state. Model 6 uses BB-I-BPL. Re- 
placement of BPL with CompTT does not change our 
conclusions for the properties of the BB component and 
the fraction of Comptonization as long as the cold-seed- 
photon solution is used. 

There is a question of whether we could also see a ther- 
mal disk in the hard state, where the trend from Figure[5] 
would imply fcTmcd ^ 0.5 keV. Such disks have in fact 
been observed in the hard state o f a number of BH X - 
ray transients ([Miller et al.ll2006al lbl; [Rvkoff et al.ll2007D . 
with instruments that have better sensitivity at low en- 
ergies than RXTE. However, being able to detect these 
components in NSs is not only a matter of low-energy 
sensitivity. One also needs a relatively low interstellar 
absorption. Moreover, in NSs a disk would not be the 
only thermal component contributing around 0.5 keV, 
since there is also emission from the boundary layer, and 
it is difficult to disentangle the contribution from the two 
components at low luminosities when the spectrum is not 
dominated by these thermal components. 

If the boundary layer is an equatorial belt of BB 
emission, then Model 6 implies that the visible sur- 
face area of the belt remains nearly constant, with 
Nhh/Nhurst ^ 1/16, over a wide range in Lx (0.005- 
0.5 -ZjEdd)- Current theories predict that the geometry 
and physical processes in th e boundary layer sh ould vary 
significantly with Lx (e.g . . Ilnogamov fc Su"iiv aev 199!|; 
iPjppham fc SunvaevI [200l iKluzniak fc Wagonen il9850 . 
and some obser vations have been interpreted in such a 
manner (Chur ch et al.|[200 2[). Our results are much sim- 
pler than expected, suggesting that a geometrically thin 
accretion stream impacts a rather well-defined surface 
area, where the gas radiates efficiently before spreading 
over the NS. This scenario would seem to require that the 
NS lies within its ISCO, whic h pinches the vertical str uc- 
ture of the accretion stream ([Abramowicz et al1ll978[ ) as 
it flows toward the NS. This topic must be investigated 
in further detail. 
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TABLE 1 

X-RAY Sources and Observations Prior to 2006 January 1 





number of total/used 


Time of total/used 


Nil 


Distance 


Source name 


observations 


observations (ks) 


(10^2 cm-2) 


(kpc) 


Aql X-1 


393/333 


1402/1252 


0.5 


5 


4U 1608-52 


523/459 


1310/1204 


1.0 


3.6 



TABLE 2 
The spectral models 



Hard State Soft/Transitional State 



Model 


Description 




Problems'^ 


Description 


xU-)'" 


Problems'^ 


1 


BPL-fMCD 


0.93(0.21),0.88(0.16) 


R 


BPL+MCD 


0.99(0.38), 1.07(0. 36) 


R L P 


2cold 


CompTT-l-MCD 


1.00(0.28),0.90(0.17) 


R 


CompTT-l-MCD 


1.04(0. 39), 0.93(0.25) 


R L P 


2hot 


CompTT-l-MCD 




R 


CompTT-fMCD 




L P 


3 


BPL-I-BB 


0.94(0.21),0.90(0.17) 




BPL-I-BB 


1.06(0.41), 1.24(0.49) 


P 


4cold 


CompTT+BB 


0.99(0.24),0.93(0.19) 




CompTT-l-BB 


0.95(0. 33), 0.90(0.25) 


P 


4hot 


CompTT-l-BB 




T 


CompTT-l-BB 




T P 


5d 








MCD-I-BB 


1.23(0.36), 1.17(0.38) 




6 


BPL-I-BB 


0.94(0.21),0.90(0.17) 




MCD+BB-I-CBPL 


1.15(0.30), 1.02(0.27) 





Note. — All models also include an interstellar absorption component and a Gaussian line. The notations BPL, MCD, BB and 
CompTT refer the bknpower, diskbb, bbodyrad, and comptt models in XSPEC, respectively. CBPL is a constrained BPL with the break 
energy E^^ fixed at 20 keV and the initial photon index Fi forced to be smaller than 2.5. means that the information is not available 

or meaningless because the corresponding model does not work. " — " means that none of the four problems listed in ^4.2.2| is found to 
apply to the corresponding model/spectral state. 

^ The subscripts cold and hot denote the cold-seed-photon models and hot-seed-photon models, respectively.^ The mean xt standard 
deviation. The two columns are for Aql X-1 and 4U 1608—52, respectively.*^ See ^4. 2. 2 1 for the meanings of these characters.'^ Only 
observations in the soft state are fitted with this model. The mean xt for this model is from observations with source intensity > 800 
counts/s/PCU. 



TABLE 3 

BEST-FITTING PARAMETERS OF TWO SAMPLE SPECTRA 



Data Model kTs keV kT^ keV r AjT^cd^^^bb keV A^mcd,A^bb x5(d.o.f) 



Hard-state sample 


CompTT-l-MCD(cold) 


< 0.35 




-0.36 


1 70+0.07 


1 7+0.3 
^•'-0.3 


1.49(78) 


Hard-state sample 


CompTT-l-MCD(hot) 


^■J^"-0.12 


i5.5ll:g 


K ^4+0.32 


1 00+0 05 


44.8i;s:? 


1.32(78) 


Hard-state sample 


CompTT-l-BB(cold) 


< 0.4 


i5.olJ:° 


^■"^•-'-0.14 


, r,i+0.04 

^•^^-0.04 




1.25(78) 


Hard-state sample 


CompTT-l-BB(liot) 


^■^•^-0.18 


16.6+11 


2 24+°-" 


n 67+0 05 
^•"'-0.05 


227.9t^t;J 


1.59(78) 


Soft-state sample 


CompTT-l-MCD(cold) 


< 0.5 


2 5+°-l 


16.6+2'5 




68.7t5^-j 


1.27(53) 


Soft-state sample 


CompTT-l-MCD(hot) 


^■^^-0.03 


2 6+°-l 
^•"-0.1 


10 1+0-3 


<0.7 




0.95(53) 


Soft-state sample 


CompTT-l-BB(cold) 


< 0.5 


2 5+°-^ 


6.08t0;g 


1 24+0-04 


164 ^+^'^■^ 

i04.i_2g 4 


1.07(53) 


Soft-state sample 


CompTT-l-BB(hot) 


1 n2+°03 

^■""^-0.04 


2 6+°-l 
^•"-0.1 


^■""-0.16 


£ 0.5 




0.95(53) 



Note. — The detailed information of these two spectra are given at the beginning of !|4.2. 1| and their unfolded spectra are 
shown in Figures 4 and 5. The notations "cold" and "hot" refer to the cold and hot-sced-photon models respectively (see ^4.2. H . 
For CompTT-|-MCD models, a spherical geometry was assumed for CompTT component, while for CompTT-fBB models, a disk 
geometry was assumed instead. N-^xcd. and A^^bb are normalizations (in unit of km^ at distance 10 kpc) of MCD and BB, respectively. 



